Entanglement in atomic resonance fluorescence 
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The resonance fluorescence from regular atomic systems is shown to represent a continuous source of non- 
Gaussian entangled radiation propagating in two different directions. For a single atom entanglement occurs 
under the same conditions as squeezing. For more atoms, the entanglement can be more robust against dephasing 
than squeezing, hence providing a useful continuous source for various applications of entangled radiation. 
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The resonance fluorescence of a single atom played an 
outstanding role when searching for a radiation source that 
clearly displays the quantum nature of light. A single two- 
level atom, after being excited by a driving laser, can only emit 
a single photon at once and then it must be re-excited, before 
another photon can be emitted. This feature led to the pre- 
diction of the nonclassical effect of photon antibunching [ 1 1, 
which could first be demonstrated in the resonance fluores- 
cence of an atomic beam [2 1 and later with a trapped ion Q. 

Sub-Poissonian photon statistics could also be observed for 
the first time in resonance fluorescence [4J and squeezing in 
resonance fluorescence was predicted |5|. Even the resonance 
fluorescence of many atoms can show squeezing, which re- 
quires stable phase relations between the emitters. This can 
be achieved by a regular arrangement of the atoms |6|, or by 
detection of the fluorescence in the forward direction with re- 
spect to the pump beam Q. Squeezing from strongly driven 
regular atomic systems has also been studied [8|. In an early 
experiment with regular atoms, the interference of the fluores- 
cence of two trapped ions was demonstrated [9 |. Squeezing 
in resonance fluorescence could be observed for samples of 
many atoms iflOl . 

A direct demonstration of squeezing in the resonance fluo- 
rescence of a single atom has not been realized yet. For this 
purpose it was proposed to apply homodyne correlation mea- 
surements 1 11]. The method has been further developed to 
detect general field correlation functions by balanced homo- 
dyne correlation techniques II 121 . This opens possibilities to 
study the most general nonclassical features of the atomic res- 
onance fluorescence radiation [13]. As an example, intensity- 
field correlation functions could be measured in resonance flu- 
orescence 1 14 1, which is a first step in such a direction. 

In the context of applications for quantum information pro- 
cessing, among the manyfold of nonclassical effects entan- 
glement became of particular importance, for recent reviews 
see lfl"5l[l6l . A variety of possibilities to create entanglement 
in atoms, e.g. via cooperative fluorescence, has been stud- 
ied, see ifTTll and references therein. A cold atom in a cavity 
can serve as a stable source for entangled EPR-type photons, 
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by utilizing the coupling between the scattered photons and 
the quantized atomic center-of-mass motion [ 18 1. To our best 
knowledge, however, the resonance fluorescence radiation it- 
self has not been studied yet with respect to entanglement. 

In the present Letter we study the resonance fluorescence of 
a regular system of atoms, separated from each other by sev- 
eral wavelengths of the driving laser. Entanglement is found in 
the electric field strengths of the fluorescence radiation prop- 
agating in different directions from the atomic sample. Re- 
markably, the robustness of the entanglement against dephas- 
ing of the atomic transitions increases with increasing number 
of atoms. This renders it possible to design continuously ra- 
diating sources of non-Gaussian entangled light for various 
applications. 

For simplicity, we deal with the resonance fluorescence 
of atomic two-level systems. The theoretical background is 
given in several textbooks, e.g. fl9l . Let us consider an en- 
semble of N non-interacting two-level atoms, located at the 
positions (n = 1 . . . N). The electric field strength oper- 
ator at position r and time t is given in source-quantity repre- 
sentation as 

JV 

E( ± )(r,t)=E| ±) (r,t)+^EW(r,i) ; (1) 

n=l 

where the indices 'f and V denote the free and source 
field parts, respectively. The superscript (±) indicates posi- 
tive/negative frequency field components. For our purposes, 
the geometry is chosen such that the pump laser does not hit 
the detectors, hence the corresponding free-field contribution 
is in the vacuum state, (• • • e[ + ') = (e| - ' • • • )=0. In the nor- 
mally ordered correlation functions considered in the follow- 
ing, we can throughout omit the free field contributions in 
Eq. ([T]i and only keep track of the source fields. 

The source field part of the resonance fluorescence field of 
a two-level atom in rotating wave approximation is given by 

E(+)(r,t)=g(r-rW)i^ ) (i(")) ) (2) 

with A$ = \a^)(b^\ ({a,6}=l,2) being the atomic 
flip operator of the n-th atom at the retarded time t^ n ' — 
t— |r— r^|/c. The function g(r— r^"-') relates the atomic 
source operators to the source field parts of the radiation field. 
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Assuming the scattering center of the atoms at r = and ap- 
plying the far field approximation, g reads as 



v g(r) = 



d (d • r)r 



Idllr 



2 " 



(3) 



where d is the atomic dipole moment and v a scaling factor. 

The atomic flip operators obey the equal time rela- 
tion A < ^A^=Sb c A^, with 5 being the Kronecker-5. Their 
quantum averages are obtained from the optical Bloch equa- 
tions for the density matrix of a single atom. One has to 
include phase factors exp(—i<Pj ), according to (6), to ac- 
count for the interference between all atoms and the positions 
of the detectors. Here n (n = 1, . . . , N) numbers the atoms 
and j labels the positions rj of the used points of observa- 
tion. Since the atoms are separated by linear distances large 
compared to the pump-laser wavelength, the correlation func- 
tions of products of atomic operators of different atoms fac- 
torize, (A^A^)={A^){A^) for m^n. That is, under 
such conditions cooperative effects are negligible. 

In the balanced homodyne correlation measurements under 
consideration, the interference of the signal field with the lo- 
cal oscillator is the observed signal, for details see fl2l . The 
recorded correlation functions contain the slowly-varying sig- 
nal field operators, (r, t) = (r, t)e ±J " lot , where uj 1o 
is the frequency of the local oscillator. The fields at the points 
of observation r,- (in the following we choose j = 1,2) can 
be related to the Bosonic operators a,j t \ as 



E(+)(r j)i ) = ^f A (r J )a i , A (t). 



(4) 



The mode functions describe those fields which are 

recorded by homodyne correlation measurements within the 
spatio-temporal resolution of the detectors. 

An important class of entangled states can be characterized 
by the negativity of the partial transposition of the quantum 
state under study [20|. The Peres condition for entanglement 
has been reformulated as [21 1 
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in terms of the partial transposition of an operator P f, with / 
being a function of the Bosonic operators describing the sys- 
tem f2D . for the multi-partite extension see |22|. Expanding 
the operator / into a Taylor series of hj, the entanglement con- 
ditions lead to a hierarchy of negativity conditions in terms of 
minors whose entries are moments of aj and at . For a bipar- 
tite system consisting of two modes, a simple example of such 
an entanglement condition is 
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This condition can be generalized for the multi-mode fields 
defined in Eq. Now the operator / in the condition {5]) is 



chosen as a Taylor series in terms of the fields in two space- 
time points (i — 1. 2), (r^ , t), for which we use the short 
hand Ef\ We can again formulate the entanglement condi- 
tions in terms of a hierarchy of minors. The quantum-field 
theoretical counterpart of the two-mode entanglement condi- 
tion (|6) reads, in terms of the minor fi, as 
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The occurring correlations of the fields in two space-time 
points relative to the radiating source describe the measured 
behavior for the often used assumption of the detector areas 
being small compared to the spatial coherence area of the 
fields and the time resolution being smaller than the coher- 
ence time. For practical applications the corresponding fields 
can also be collimated by apt optics and propagated to two 
different regions, e.g. to Alice and Bob, which are far apart 
from each other and far from the source. Both these possibili- 
ties, the detection and the propagation of the entangled fields, 
are specified by choosing the mode functions in Eq. Q to de- 
scribe localized and propagating fields, respectively. 

Expanding the minor fi in Eq. Q, one readily gets as con- 
dition for entanglement 



(AE^ ) AEi +; )(AEf J AE^ ; )-|(AEf ; AE. 

Here we have introduced the symbol AO = O — (O) for an 
arbitrary operator O. If one would further reduce the minor 
in the condition (|7]i to the dimension 2x2, it is easy to see 
that the resulting minors are always non-negative. Hence the 
entanglement condition (jH} under study is a rather simple ex- 
ample, but it will turn out in the following to be sufficient for 
demonstrating entanglement in resonance fluorescence. 

For the resonance fluorescence from identical two-level 
atoms, the dependence on the number N of atoms in // is in 

(n) 

the phase factors tp\ . After some straightforward algebra, 
the left hand side of the condition ([8) reads as 
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where Cba={A a b) are the slowly varying density matrix ele- 
ments of the identical two-level atoms and 1? is the angle be- 
tween the two functions gi and g 2 , with gj = g(ij-). For 
random positions of the atoms, the phase combinations in the 
second and third line in Eq. |9]l are averaging to zero, so that 
only the first line contributes. Under these conditions the mi- 
nor is always positive semidefinite, and hence entanglement 
does not exist or it cannot be identified by the condition d8). 
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Let us consider now the situation for a regularly arranged 
sample of N atoms, whose resonance fluorescence is known 
to show squeezing J6), contrary to the situation for randomly 
distributed atoms. Under optimal conditions all the phase dif- 
ferences, p!? 1 ' —ipj , vanish. This requires that all atoms lie 
on the central plane between the two detector positions, and 
the (not necessarily equal) distances between the atoms are 
integer multiples of the pump laser wavelength. In such an 
optimal configuration, the condition ([8]) reduces to 



- 2 ~ 2 sin 2 (i?) 
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It is straightforward to reduce the condition (jTOj for the gen- 
eration of entangled spatial light modes by TV atoms, by using 
the general property O22>|cr2i| 2 >0, to 
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Inserting the stationary values of <jy from |fl9l , we rewrite the 
entanglement condition in terms of the atomic and driving- 
field parameters and obtain for the Rabi frequency £!r 
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where T% and r 2 are the energy and phase relaxation rates, 
respectively. Together with the deviation A of the laser fre- 
quency from the atomic transition, these rates determine the 
maximal driving field for which entanglement is obtained. 

The Rabi frequency is chosen real, so that the right hand 
side of the inequality ( 12 1 has to be positive. For N=l, this 



condition is equal to the requirement for the possibility to ob- 
serve squeezing in single atom resonance fluorescence J3]. 
For iV>l the driving field can be increased by a factor of 
about y/~N compared with a single atom. More importantly, 
with an increasing number of atoms the ratio T2/T1 of the 
individual atoms can exceed the value of one, being the ulti- 
mate limit beyond which squeezing in resonance fluorescence 
of a regular TV-atom system disappears [6|. Thus, for such 
a system entanglement is more robust against dephasing than 
squeezing. 

We would like to emphasize at this point that the reso- 
nance fluorescence not only provides a continuously radiat- 
ing entangled source with high stability against dephasing, 
but it is also non-Gaussian. This is easily verified, since cor- 
relations containing operator products m are zero 
when n+m exceeds the number N of atoms. It is important 
that non-Gaussianity is "at times necessary to successfully 
perform quantum information tasks", cf. [23 ] and Refs. [4-17] 
therein. This includes applications for entanglement distilla- 
tion and swapping, quantum teleportation, universal quantum 
computing, and others. That is, for some important quantum 
tasks the non-Gaussian property is a desired resource. For 
this reason, protocols have been established for the aim to 
de-Gaussify intially Gaussian entangled radiation fields, see 
e-g- EH. The resonance fluorescence source under study is 



per se non-Gaussian and hence does not require such addi- 
tional procedures. This opens a wide field of possible appli- 
cations whose study is beyond the scope of this Letter. 

In Fig. uj we show the minor fi normalized by ^ _4 A^ 2 cr| 2 , 
cf. Eq. ^ together with Q, for a linear chain of atoms. 
All atoms are exactly aligned along the z-axis, with an inter- 
atomic distance of ten wavelengths of the driving laser. The 
direction of the laser is chosen to be almost parallel to the z- 
axis, so that the Rabi frequencies are kept nearly constant for 
a sufficiently large number of atoms. The two detectors are 
positioned in the far field region (at distance r from the coor- 
dinate center), detector 1 in the y-direction and detector 2 can 
be moved in the rr-y-plane. Hence, the position of the latter 
only depends on the azimuth angle fa. The system parameters 
are chosen such that the incoherent part of the atomic fluores- 
cence plays a significant role, 022/|cr2i| 2 =3.5. We remind the 
reader that under these conditions squeezing in resonance flu- 
orescence would not occur. According to the condition (Hi, 
entangled light is irradiated by three or more atoms. We con- 
sider the situation for the numbers of atoms being N—2, 3, 4. 
It is clearly seen from Fig. [T] that entanglement occurs for 
N > 3. Even for the chosen relatively poor coherence prop- 
erties, entanglement of the radiation fields in two directions is 
found in almost every direction. Note that for opposite direc- 
tions of the detectors relative to the atomic sample (ri=— r 2 , 
i.e. sin?9=0) the entanglement vanishes in general. 





0.05 











b 

(S 


-0.05 


1 


-0.1 


a. 






-0.15 




-0.2 




FIG. 1: The normalized minor /i is shown, as a function of the angle 
<f>2, for a linear chain of atoms with N = 2 (solid), 3 (dashed), and 4 
(dotted). The pumping strength yields o"22/|o"2i| 2 =3.5. 

The linear chain cannot be experimentally realized to arbi- 
trary precision. For more and more atoms the negativity of 
the minor increases together with the robustness of the corre- 
sponding entanglement against dephasing allowing less accu- 
racy without destroying the entanglement. For the purpose of 
describing a realistic system, let us consider a linear Paul trap 
for a chain of ions. Note that present technologies, such as the 
multizone trap-array architectures designed for quantum com- 
putation [25 1, may offer alternative possibilities to realize an 
entangled radiation source of the type under study. The posi- 
tions of the ions are rather well fixed in the linear trap. How- 
ever, they are not perfectly equidistant for more than three 
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ions, due to the distance dependence of the Coulomb interac- 
tion within the chain of ions, for more details see e.g. Il26l . 
The atomic distances are scaled by an externally controlled 
length parameter 7, which can be properly adjusted to the 
number N of trapped ions. 

One also cannot ignore fluctuations in the positions, as the 
individual ions are not really fixed due to the uncertainty prin- 
ciple. Assuming a harmonic potential around each of the equi- 
librium positions of the ions, the frequency of these harmonic 
oscillators are nearly equal to the frequency of the trap oscil- 
lation. This simple assumption overestimates the variances of 
the atomic positions, since the potential around a given ion 
rises quickly due to anharmonic effects. Hence the true posi- 
tion uncertainty can be much smaller than that of an ion in the 
full trap potential. Nevertheless, we choose this larger value, 
to effectively account for other possible experimental imper- 
fections. Using elementary quantum mechanics and statistics, 
we can approximate the physical limit of the variance of the 
individual ions as 



Az 
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uous radiation source emitting two multi-mode light beams 
in different spatial directions, which show bipartite entangle- 
ment. The desired radiation source consists of the atomic 
resonance fluorescence of a regular system of several non- 
interacting atoms. This represents a non-Gaussian source of 
entangled radiation. For a single atom, entanglement is ob- 
tained under the same conditions as required for the realiza- 
tion of squeezing in resonance fluorescence. For a system of 
N atoms, the situation can be significantly improved, as by in- 
creasing the number of atoms the entanglement of the emitted 
light beams is getting more robust against atomic dephasing. 
As an example, we have estimated the practical limitations of 
the achievable entanglement for a system of ions in a linear 
trap. Altogether, the resonance fluorescence of regular atomic 
systems is a continuous and non-Gaussian source of entangled 
light, which may open interesting perspectives for various ap- 
plications in quantum information technology. 

This work was supported by the Deutsche Forschungsge- 
meinschaft, SFB 652. 



(13) 



where M is the mass of one of the ions and Q its charge. As- 
suming for example Q = e and M = 3.3309 x 10~ 25 kg 
and an optical transition wavelength A = 194.2 nm, we ob- 
tain Az 



.3/4 



This is the situation to 



\z w 1.014 x 10- y m|_ 
be realized with mercury ions, cf. the experiments in (9)- As 
the variance should be small compared to the transition wave- 
length to keep statistical averaging effects in the minor /1 given 
by Eq. (|9]l sufficiently small, we assume a limit of Az<0.1A. 
This leads to a maximal value for 7 of about 50A. Thus we 
conclude that linear traps would be suited to provide sufficient 
positioning accuracy of the atoms to realize a continuously ra- 
diating entangled-light source with regularly arranged atoms 
for average distances large compared to the wavelength. 

To verify these theoretical predictions by simulations, we 
choose the parameter 7 in such a way, that for each N the mi- 
nor fi becomes minimal. The resulting Az from Eq. ( 13 1 is 



implemented as an upper bound for a random deviation from 
the ideal position along the z-axis, variations perpendicular 
to this axis are assumed to be negligible due to properly cho- 
sen potentials. The minors have been calculated for the same 
number of ions and the same parameters as in the ideal case 
shown in Fig.[T] The calculations for the realistic chain are in 
reasonable agreement with the idealized results given above. 
For N = 2,3 the uncertainties of the ion's positions mod- 
ify the idealized result by less than one percent. For four 
ions in the linear trap, we could optimize the negativities to 
about 98.8% of the ideal case. Thus it is possible to obtain 
almost optimal entanglement in the atomic resonance fluores- 
cence with ions in a linear trap, as long as the number N of 
ions is not too large. By applying alternative technologies for 
realizing regular atomic systems, e.g. by trap arrays EBI . one 
may also overcome this limitation. 

In conclusion, we have studied the realization of a contin- 
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